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We study the formation of domain walls in a phase transition in which an 5*5 x Z2 symmetry 
is spontaneously broken to 5*3 x 5*2. In one compact spatial dimension we observe the formation 
of a stable domain wall lattice. In two spatial dimensions we find that the walls form a network 
with junctions, there being six walls to every junction. The network of domain walls evolves so 
that junctions annihilate anti-junctions. The final state of the evolution depends on the relative 
dimensions of the simulation domain. In particular we never observe the formation of a stable 
lattice of domain walls for the case of a square domain but we do observe a lattice if one dimension 
is somewhat smaller than the other. During the evolution, the total wall length in the network 
decays with time as t""'^^, as opposed to the usual t~^ scaling typical of regular Z2 networks. 

PACS numbers: 98.80.Cq, 05.70.Fh 



I. INTRODUCTION 

Topological defects are routinely observed in con- 
densed matter systems and may be expected to form 
during phase transitions in the early universe. Observa- 
tional constraints on various types of defects have already 
played a major role in the development of cosmology - the 
idea of inflation was introduced in large part to solve the 
monopole problem. Recently, much attention has been 
devoted to cosmological models motivated by String/M 
theory. Symmetry breaking patterns in these models can 
be quite complex with the production of topological de- 
fects being a generic phenomenon. 

It has been known for some time that, while topologi- 
cal considerations are sufficient for proving the existence 
of defects, the properties and interactions of defects de- 
pend on the details of the particular model. These issues 
were highlighted in earlier studies of domain walls in an 
SU{5) X Z2 model. In contrast to the commonly stud- 
ied Xcj)'^ model, the vacuum manifold now consists of two 
disconnected pieces, each of which is a 12 dimensional 
continuum. For topological reasons, kink solutions must 
exist. However, a 12 dimensional continuum of possibili- 
ties exist when determining the boundary conditions that 
the minimum energy kink solution satisfies. Topological 
arguments alone are not sufficient in order to determine 
these boundary conditions and a more elaborate analysis 
is needed if one wants to find the kink solution ^ . Only 
once the boundary conditions are known, can the explicit 
solution be constructed and the properties worked out. 
In the case of domain walls in SU (TV) x Z2 models, this 



exercise has already yielded some surprises. For example, 
it was found that the symmetry group inside the core of 
a domain wall is generally smaller than that of the vac- 
uum P, 13, 0, Q . The interaction of kinks in these models 
shows that kinks and antikinks can repel as well as at- 
tract ^ . In Ref . it was argued that these features are 
in fact generic when large group symmetries are involved. 
Solutions of similar type were also discussed in . 

The formation of kinks and domain walls in phase tran- 
sitions has been studied in earlier work but most of this 
work only deals with the simplest of systems, such as the 
Z2 kink in a Xcj)'^ model. Based on the expected initial 
density and scaling of these walls, cosmological impli- 
cations have been derived. However, as we show here, 
the properties of a domain wall network in more complex 
models can be dramatically different from that of the wall 
network in the Xip^ model. Earlier work along these lines 
can be found in Ref. ^ where they consider the forma- 
tion and evolution of walls, and in Refs. '5,1^ where 
the authors study the fate of walls in 0{N) motivated 
models. 

An example of new physics that one can expect when 
considering domain wall formation in these more complex 
models was provided in . There it was shown that the 
SU{N) X Z2 models allow for the existence of domain 
wall lattice solutions. The key ingredient is the repulsion 
between kinks and antikinks. The lattice is a periodic 
sequence of repelling walls and antiwalls that are parallel 
to each other. It is more challenging to find a model 
in which the lattice is stable. In 11], a model in which 
S5 X Z2 {Sn is the permutation group of n objects) breaks 
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to 5*3 X S2 was used as an example of a model allowing 
for stable domain wall lattices. It was argued that in 
one spatial dimension, in the limit of many correlation 
domains, the probability of forming a lattice tends to 
unity. 

In this paper we follow on the work in and nu- 
merically investigate the formation of lattices during re- 
alistic phase transitions in (1-1-1) and (2-1-1) dimensions. 
Our results in 1-1-1 dimensions corroborate the argument 
in and kink lattices are observed to form with near 
certainty. The simulations in 2-|-l dimensions, however, 
yield an unexpected feature that, with hindsight, might 
have been anticipated from the discussion in Ref . . In- 
stead of forming a lattice, the domain walls form a net- 
work with junctions. Six walls meet at a junction and the 
relaxation of the network is controlled by the dynamics 
of the junctions. 

We find that this leads to considerable changes in 
the long time dynamical properties of the wall network. 
Whereas the total length of a Z2 network of domain walls 
is expected to decay as , our junction dominated net- 
work has a lower scaling power: t^^''^^ . Junction motion 
and junction/anti-junction annihilation processes clearly 
slowdown the long time evolution of the network. We did 
not observe the formation of a lattice of domain walls as 
the final evolution state in our simulations on a square 
spatial grid (with periodic boundary conditions). How- 
ever, when the spatial grid is rectangular, with one di- 
mension somewhat smaller (by roughly a factor of 3) than 
the other, we did observe kink lattice formation, even in 
2-1-1 dimensions. 

The paper is organized as follows. In Section ^1 we 
introduce the model and review its features relevant for 
the process of domain wall lattice formation. Section 
mil contains the simulation results in both one and two 
spatial dimensions, including a discussion of the scaling 
regime of the wall network in the 2-1-1 case. Finally 
we discuss how these results may vary for more general 
systems in different dimensions. 



II. REVIEW OF DOMAIN WALL LATTICES 

In this section we review some of the results from 
Refs. 0, 0, 0; IE 0| relevant to domain wall lattices. 
In ^3 Sr, X Z2 was used as an example of a model in 
which such lattice solutions are stable. However, the do- 
main wall solutions themselves were identical to those in 
SU{b) X Z2 and we prefer, for clarity reasons, to use this 
model to describe their main properties ^. 

Consider an SU{5) x Z2 field theory described by a 
Lagrangian: 



L = Tr(9^$)2 - y($) 



(1) 



^ Most of the results in Sectionmcan be generalized to SU(N) X Z2 
with N > 3. The interested reader is referred to 



where $ is an 5*17(5) adjoint and V{^) is invariant un- 
der SU{5) X Z2- Let V{^) be such that the expectation 
value of $ spontaneously breaks the symmetry down to 
SU{3) X SU{2) X U{l)/Z3 X Z2. We will choose V{<P) to 
be a quartic polynomial: 

T/($) = ~m^TT[<P^] + /i(Tr[$2])2 + ATr[$'*] + Vq (2) 

where Vo is a constant chosen so that the minimum of 
the potential has V — 0. The Lagrangian is symmetric 
under $ ^ — $ and it is the breaking of this Z2 symmetry 
that gives rise to topological domain wall solutions. The 
desired symmetry breaking is achieved in the parameter 
range 



h _T_ 



(3) 



The vacuum expectation value (VEV), <I>o is (up to any 
gauge rotation) 



$0 



with r\ — mj ^/JJ and 



A' 



diag(2,2,2,-3,-3) 



(4) 



(5) 



In Refs. 0,1^1^ it was found that there are several kink 
solutions in this model corresponding to different choices 
of asymptotic field configurations. The necessary condi- 
tion for the existence of a kink solution $^(2;), proved in 
Ref. Q, is [$fc(x = ±oo),$/c(x)] = 0. That is, the so- 
lution must commute with its asymptotic values. It was 
also proved in Ref. l] that, when searching for SU{5) 
kink solutions, one can work in the Cartan subalgebra 
of SU{5), which is equivalent to restricting $fe(x) to a 
diagonal matrix form: 

$(x) = /i(a;)A3 + f2{x)Xs + f3{x)T3 + U{x)Y , (6) 

Ag, T3 and Y are the diagonal generators of 



where A3 
SU{5): 



A3 
As 

T3 
Y 



idiag(l,-l,0,0,0) , 
diag(l, 1,-2, 0,0), 



2V3 



-diag(0, 0,0, 1,-1) 



1 



2V15 



diag(2,2,2,-3,-3) 



(7) 



As shown in Refs. Q, |^ |^, the kink solution with 
least energy is achieved if (up to global gauge rotations) 
^{~oo) = $_ = $0 and 



$(-1-00) 



$4 



diag(2,-3,-3,2,2) (8) 
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The minus sign in front of in Eq. ^ puts $+ and 
in disconnected parts of the vacuum manifold. Also, 
two blocks of entries of <!>+ are permuted with respect to 
those of In other words, <&_ and — are related 
by a non-trivial gauge rotation. The kink solution (or, 
domain wall solution, in more than one dimension) can 
be written down explicitly in the case when h/X = —3/20 
[2, 3]: 



1 — tanh((Tx) 



1 + tanh(cra;) 



$4 



(9) 
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where a ~ m/ \f2. For other values of the coupling con- 
stants, the solution has been found numerically ^J. 
The topological charge of a kink can be defined as 



(10) 



7? 



where <J>fl and are the asymptotic values of the Higgs 
field to the right (i?) and left (V) of the kink. (The rescal- 
ing has been done for convenience.) Then the charge of 
the kink in Eq. Q is: 



Q« =diag(-4,l,l,l,l) 



(11) 



Similarly, one can construct kinks with charge matrices 
Q(*^ {i = 1,...,5) which have —4 as the ii entry and 
+1 in the remaining diagonal entries. Hence there are 
kink solutions with 5 different topological charge matri- 
ces. Individually, the kinks can be gauge rotated into 
one another. But when two kinks are present, the differ- 
ent charges are physically relevant. This is most easily 
seen by noting that the interaction between a kink with 
charge Q^*^ and an antikink with charge Q'--'^ = — Q*--'-' is 
proportional to Tr(Q(*)g(j)) 0. Then we have 

Tr(gWg(^)) = -20 ff i = J 

= +5ifz^i (12) 

The sign of the trace tells us if the force between the kink 
and antikink is attractive (minus) or repulsive (plus). 
Hence the force between a kink and an antikink with dif- 
ferent orientations (i ^ j) is repulsive. This observation 
is key to the construction of kink lattices. 

A kink lattice is a periodic sequence of kinks with 
charges such that the nearest neighbour interactions are 
repulsive. One can write down a sequence of charges that 
can form a kink lattice 0|: 

...Q(i)g(^)g(3)Q(i)Q(^)Q(3^.. (13) 

and the sequence just repeats itself. This sequence is the 
minimum sequence for which the nearest neighbour in- 
teractions are repulsive. Another way to write the kink 
sequence is to write it as a sequence of Higgs field ex- 
pectation values. We write this sequence for the above 
minimal lattice: 

... ^ +(2, 2, 2, -3, -3) ^-(2, -3, -3, 2, 2) 
^ +(-3, 2, 2, -3, 2) ^-(2, -3, 2, 2, -3) 
^ +(2,2,-3,-3,2) ^ -(-3,-3,2,2,2) 
^ -K(2,2,2,-3,-3)->... (14) 



The minimal lattice of 6 kinks is not the only possibility. 
A sequence of 10 kinks in the N = 5 case is aesthetic in 
the sense that it uses all the 5 different charge matrices 
democratically: 

There can be longer sequences as well. 

A detailed stability analysis in Ref. revealed that 
all domain wall lattices in SU{5) x Z2 are unstable. 
For example, the lattice in Eq. H13|l has three unsta- 
ble modes, corresponding to rotations in the 1-3, 1-5, 3-5 
blocks. The instability comes from the fact that an iso- 
lated kink has zero modes corresponding to rotations in 
field space - for example, a kink with charge Qi can be 
rotated into the kink with charge Q3 without any cost 
in energy. When a kink of charge gi is placed near an 
antikink of charge Q3, the zero mode becomes an unsta- 
ble mode, making it favourable for gi to rotate into Q3 
after which the kink and antikink can annihilate. 

To illustrate that stable lattices can exist, one can sim- 
ply start with the model in which the zero modes are com- 
pletely absent right from the start. As in Ref. ^3|> us 
consider the model of four real scalar fields fi {i ~ 1, .., 4), 
with 



1 4 



(16) 



and 



V = — — 



4 

E 

i=l 



f! + 



4 



1=1 



3 



ffr + i} fa 



1 



7_ 

30 



fi 



(17) 



The fields /,; are defined as in Eq. lO and this model 
has been obtained by substituting Eq. © into Eq. 1^. 
This four field model does not have the continuous SU (5) 
symmetry of the model in Eq. The only remnant of 
the SU (5) symmetry corresponds to the permutation of 
the five diagonal entries of $. In addition, the model also 
has the Z2 symmetry under which fi — fi- Hence the 
model has an S5 x Z2 symmetry. 

A vacuum of the model is given by /i = = /2 = 
/a and /4 ^ 0. This breaks the symmetry to £"3 x 6*2, 
corresponding to permutations of $ in the SU{3) and 
SU{2) blocks. The vacuum manifold consists of 5! x 
2/3! X 2! = 20 discrete points. If we fix the vacua at a; = 
—00, this implies that there are 20 kink solutions in the 
model. All these 20 kink solutions have been described 
in Ref. ^. 

The construction of kink lattices proceeds exactly as 
in the SU (5) case above because the off-diagonal compo- 
nents of $ vanish there. Hence the S^ x Z2 model contains 
kink lattice solutions as well. Furthermore, these lattices 
are stable because the dangerous rotational perturbations 
are absent by the very construction of the model. 
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III. LATTICE FORMATION 

There are more repelling kink-antikink pairs in S5 than 
attracting ones. Hence, it is reasonable to expect that af- 
ter the phase transition, all attracting walls will eventu- 
ally annihilate and the remaining walls will be repelling. 

The probability of forming a domain wall lattice in 
one spatial dimension for the model was estimated 
in to be unity if the total number of kinks formed 
is large. Numerical simulations, presented below confirm 
this expectation. In two spatial dimensions with periodic 
boundary conditions, our numerical results do not show 
lattice formation. Instead the walls form a network with 
junctions that gradually dilutes due to the annihilation 
of junctions. 

A. Numerical implementation 

In order to generate feasible initial conditions that 
may lead to the formation of domain wall lattices, we 
use a Langevin type equation based on the Lagrangian 
Eq. (|16|l . Each field will be propagated according to its 
usual equation of motion with additional dissipative and 
stochastic terms added: 

(92 _ /, + d,v + D dth - r„ (18) 

D is the dissipation constant and the stochastic force 
Ti{x,t) is a Gaussian distributed field characterised by a 
temperature T: 

2D 

{V,{x,t)V,{x',t')) = —5,.^5{x-x')5{t-t') (19) 

The amplitude of the noise in Eq. H19() is chosen so as to 
guarantee that independently of the initial field configu- 
ration and of the particular value of the dissipation, the 
system will always equilibrate towards a thermal distri- 
bution with temperature T. 

In both the ID and 2D cases, we will start by evolving 
Eq. H18|) until thermal equilibrium is reached. At that 
point we quench the system to zero temperature by set- 
ting the stochastic term, F^, in Eq. p8|l to zero. The 
fields will then settle towards the minima of the poten- 
tial and a network of domain walls will form, separating 
regions where the fields were initialy uncorrelated. Note 
that the correlation length ^ at thermal equilibrium de- 
pends on T (at high temperatures, ^ typically decreases 
with T). This allows us to have a degree of control over 
the correlation length of the fields before the quench, and 
hence over the number of independent domains that will 
form in the simulation box. This is essential if we want 
to be certain to have a number of domains large enough 
to generate a stable lattice. 

The equations of motion were discretized using a stan- 
dard leapfrog method and periodic boundary conditions 
were used. The model parameters were set to m = 
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FIG. 1: Energy contour for a ID system at different stages 
of the evolution. The top plot shows the initial thermal con- 
figuration followed by a snapshot characteristic of early times 
after the quench to T = 0. At later times the kink network 
evolves to form a stable wall lattice as shown in the bottom 
plot. 

1/(2^/6), A = 1/2 and = -3/40. The lattice spacing 
was 5x — 1. and 6t — .5. Note that for the parameters 
above the wall core is resolved by more that 10 lattice 
points which should be accurate enough for the desired 
purposes. The value of the dissipation coefficient D does 
not influence the results during the stochastic stage of 
the simulation and we set it to Z) = 1.0 to ensure rapid 
thermalization. In the next two sections we will discuss 
in more detail the role of the dissipation after the quench. 

B. Simulation in (1 + 1)D 

Our goal here is to check whether a stable kink lattice 
forms in one spatial dimension as an initial "hot" config- 
uration is quenched to T = 0. Applying the procedure 
detailed in the previous section to the system, we confirm 
that this is in fact the case. In Fig. ^ we show a series of 
snapshots of the energy profile of the fields for different 
times in a typical simulation run. The first configura- 
tion corresponds to the end of the thermalization phase, 
with large amplitude fluctuations and very uncorrelated 
fields. At this point we turn the stochastic term off and 
let the system evolve for a considerably long time. In 
this case we keep the dissipation high since we are not 
interested in the details of the dynamics of the system 
but rather in observing a stable lattice as the final out- 
come. As the energy is then dissipated, a recognisable 
pattern of "proto-kinks" starts to form, as can be seen 
in the second snapshot. The next energy profile shows 
that the field has relaxed to a fully formed network of 
kinks. During further evolution, neighbouring attractive 
kink-antikink pairs annihilate, leading to a stable lattice 
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of equidistant domain walls. This final state is shown in 
the fourth plot, a succession of mutually repelling kinks 
and antikinks. The energy of these defects corresponds to 
the value predicted analytically, and direct inspection of 
the values of the four fields fi confirms that, as expected, 
we are in presence of a genuine domain wall lattice. 



C. Simulation in (2 + 1)D 

After having established that in one spatial dimension 
kink lattices can form as a consequence of a phase tran- 
sition, we will now look at the evolution of the same 
model in the (2+l)D case. Here we expect the late time 
dynamics of the system to be dominated by networks 
of one-dimensional domain walls. As a consequence of 
the symmetry content of the model, these networks will 
be considerably more complex than the ones based on 
ordinary Z2 walls. In particular, we can expect walls to 
intersect and kink-antikink repulsion to play a role in the 
evolution. 

Our first task is to find a method to identify the walls 
starting from the field values in the simulation lattice. To 
this end it is convenient to convert the fields fi{x) back 
into the original SU{5) x Z2 matrix form, $(a;). For any 
type of kink in the model, there is at least one element 
$ of the diagonal of the field matrix that changes sign 
as one crosses the defect core 0, 0, 0] . In the particular 
case of the lowest energy kinks, those described by the 
charge matrices Q*-*'' , it is easy to see that only one of the 
<l?jj's changes sign. For these cases, Tr[$'^] also vanishes 
at the defect's core. Since this does not happen for any 
of the other kink types, we have a way of distinguishing 
between the least energetic kinks and general unstable 
ones. We thus measure at each time step the total num- 
ber of zero crossings between all adjacent lattice points, 
of both the 's and of Tr[<I>'^] . For most of the evolution 
the network consists predominantly of stable kinks and 
the two quantities coincide, defining the total wall length 
in the system. 

The simulations follow the same pattern as before, 
starting with a thermalization stage after which the 
stochastic term in the equation of motion is set to zero. 
In this case, however we are interested in studying the 
long time scaling behaviour of the network. Since we are 
mostly concerned with the relativistic limit of the theory 
the dissipation term must be set to zero during the scaling 
period. Nevertheless, since the initial thermal configura- 
tion is very energetic, we keep the dissipation high for 
some time after the quench, before setting D = for the 
scaling regime. In this way, some of the excess energy 
is dissipated away, allowing the fields to relax into a well 
formed network configuration that then can start scaling. 

In all simulations, the model parameters and space and 
time discretization steps used were the same as in the one 
dimensional runs. The fields were evolved in a = 2500^ 
point grid. 

In Fig. [3 we can see a snapshot of the wall network 




FIG. 2: High density wall network for an early evolution time, 
t = 400. 
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FIG. 3: Late time network configuration t — 3200. 



for an early time in the scaling period, t = 400. The 
network is very dense with a very high number of in- 
tersections. The walls were identified using the method 
discussed above. At this time all the kinks in the net- 
work are stable and the number of zero crossings for the 
^jj's and for Tr[$'^] coincide. By direct observation of 
the field configurations, we confirmed that higher energy 
kinks do exist for earlier times which quickly decay into 
stable ones. A secondary consequence of such decay pro- 
cesses is that stable kinks forming in pairs from a single 
unstable kink tend to remain spatially correlated near 
junctions. Since the repulsive force between kinks is ex- 
ponentially small, this pairing is relatively stable and can 
be observed for considerably late times (see Fig. 

For later times most walls annihilate, loops form and 
decay into radiation and the overall density decreases. 
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Walls tend to get straighter as this is energetically more 
favourable, though the time-scale relevant for this pro- 
cess is clearly larger than the decay time of the network. 
In Fig. Owe show the network for a later time, t ~ 3200. 
The reduction in density is remarkable, with only a small 
number of intersecting walls remaining. The way these 
intersect to form nodes is quite interesting and we can see 
in particular that around each node there are always six 
incoming walls. This is related to the fact that the mini- 
mal periodic wall lattice is composed of six walls, as dis- 
cussed in Section m It is clear that for such a configura- 
tion to be stable for very long times, the walls around the 
node must repel each other. Hence the sequence of their 
charges must correspond to a stable wall lattice pattern. 
For earlier times we observe nodes with other numbers 
of incoming walls, though, as expected, never less than 
six. High index nodes are rarer because they are both 
less likely to form and more likely to cancel with other 
nodes to form a minimal six wall configuration. Even- 
tually all nodes annihilate with anti-nodes with reversed 
"vacuum orientation" , and for very long times all walls 
disappear. In all the simulations performed in a square 
lattice, the final state always had zero wall content and a 
two dimensional wall lattice never formed. Nevertheless, 
we observed that by taking one of the linear dimensions 
of the simulation domain Ly to be smaller than the other 
Lx, a regular lattice of straight parallel walls can be ob- 
tained as the outcome of the evolution. Though more 
careful simulations need to be performed, our prelimi- 
nary results suggest that as the smaller grid dimension 
increases, the final density of walls in the lattice becomes 
smaller. When Ly/L^, the ratio of the two dimensions 
becomes larger than roughly 1/3, no lattice is formed. 
This result can be better understood by imagining the 
limit when one of the torus dimensions (say Ly) goes to 
zero and the ID case is recovered. In this situation, the 
final state of the evolution should be a series of parallel 
walls crossing the torus in the j;-direction. As Ly in- 
creases, we can expected that some of these walls will be 
able to "explore" the torus in both the x and y-directions. 
Some of these will not cross the torus, getting linked to 
other walls and forming a complex network. As the net- 
work evolves, wall nodes will annihilate, some of these 
annihilations further decreasing the number of walls that 
cross the box in the y-direction. As a result, the final 
state of the system will still be a set of straight repelling 
parallel walls in the y-direction, but with a lower density. 
It is easy to imagine that as Ly increases and approaches 
Lx, the final density will decrease, up to the point where 
no lattice will form at all. The fact that the final state 
of the transition depends on the geometry of the system 
is unexpected and may have implications in cosmolog- 
ical settings in models with extra dimensions. Still, a 
better understanding of this phenomenon and the mech- 
anisms behind it is needed before further conclusions can 
be drawn. 

Another indication that the dynamics of this system 
differs fundamentaly from regular defect networks can be 



Wall Length Scaling 
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t 

FIG. 4: Wall length versus time for both the Ss x Z2 theory 
(top curve) and the single field Z2 case (bottom curve). The 
dashed curve corresponds to the total 55 x Z2 wall length 
measured by counting zeros of the diagonal elements of $(a;). 
Except for very early times this coincides with the number of 
zeros of Tr[<l>^] (solid curve). The results shown are averages 
of 25 independent realizations. In all cases, for t < 100 the 
system is evolved with both non-zero dissipation and thermal 
noise. At t = 100 we quench the system to zero temperature 
but the dissipation term is kept until t = 125 to ehminate 
excess energy. For t > 125 we set D = 0. The power-law fits 
(bold lines) to both curves were taken between t — 200 and 
t — 1350. For later times (larger than half the simulation box 
length) finite size effects become significant. 



obtained by studying its scaling properties. In Fig. ^ we 
show a log — log plot of the total wall length versus time. 
For comparison we also include a plot of the scaling of a 
regular Z2 domain wall network. This was obtained by 
evolving a scalar X(f>'^ theory in parallel with the 6*5 x Z2 
system. In both cases after the initial thermalization 
and dissipation periods, the domain wall network enters 
a scaling regime where the time evolution of the total 
wall length is well described by a power-law N ~ 
For the Z2 theory we find a = 1.06(0.06) which is in 
reasonable agreement with both previous theoretical and 
numerical predictions of a = 1. (see [T^ and references 
therein). This is to be compared with the result for the 
55 X Z2 walls. In this case we find a = 0.71(0.02), a 
considerably lower result. Clearly, the dynamics of the 
55 X Z2 network is fundamentally different from the Z2 
case. This is not surprising, taking into account that the 
process of node/anti-node annihilation must play an im- 
portant role in the dynamics of the system. That the 
overall effect is a slowdown in the scaling, could also be 
expected on the basis that the six-wall junctions feel a 
force that is predominantly isotropic. This suggests that 
systems allowing for stable nodes with higher number of 
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crossing walls and hence more isotropic, would be likely 
to display even lower scaling powers. In the limit of infi- 
nite number of incoming walls per node the system would 
eventually become static. 

IV. CONCLUSIONS 

We have numerically studied phase transitions in a 
model with ^5 x Zi symmetry breaking down to x S2 
in one and two spatial dimensions. In one dimension, 
as expected, we find that a stable domain wall lattice 
forms. In two dimension, we observed the formation of 
a complicated network of domain walls and domain wall 
junctions. At every junction six domain walls are present. 
With time, the junctions move and annihilate, and the 
network coarsens. This feature of the network dynamics 
leads to a slowdown of the scaling regime. We found that 
the total length of domain wall in the network scales as 
t^''-^^, in contrast to the ° fall off expected for Z2 
domain walls in a A^'* model. 

We have not studied the phase transition in three di- 
mensions since the problem then becomes computation- 
ally very intensive. Even then our results can be ex- 
trapolated to three dimensions, allowing us to anticipate 
certain behavior. In three dimensions we expect the wall 
junctions to be one dimensional - somewhat like strings. 
The network of walls will coarsen as the strings come 
together and annihilate. This is very reminiscent of a 
system of cosmic strings in which each string then gets 



connected to six domain walls and the evolution of the 
networks should be similar 0, 0, 0, 0| . The new fea- 
ture in the present case is the repulsion between walls 
and antiwalls and this could lead to lattice formation at 
very late times. 

Perhaps the most important application of our study 
is to the SU{5) x Z2 case as this is a realistic particle 
physics model. As we have noted in Sec. the domain 
walls in the SU (5) model are identical to those in the S'5 
model. The difference is only in their stability properties. 
Hence we might expect some of the features of the domain 
wall networks in the S'5 case to carry over to the SU (5) 
case. Exactly how much similarity will survive is hard 
to predict and may depend on model parameters. This 
remains an important problem to explore. 
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